Collective Modes in a Unitary Fermi Gas across the Superfiuid Phase Transition 
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We provide a joint theoretical and experimental investigation of the temperature dependence 
of the collective oscillations of first sound nature exhibited by a highly elongated harmonically 
trapped Fermi gas at unitarity, including the region below the critical temperature for superfluidity. 
Differently from the lowest axial breathing mode, the hydrodynamic frequencies of the higher-nodal 
excitations show a temperature dependence, which is calculated starting from Landau two-fluid 
theory and using the available experimental knowledge of the equation of state. The experimental 
results agree with high accuracy with the predictions of theory and provide the first evidence for 
the temperature dependence of the collective frequencies near the superfiuid phase transition. 
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Collective oscillations provide powerful tools to under- 
stand the physical behavior of quantum many-body sys- 
tems from different points of view and to test fundamen- 
tal theories. On one hand, collective modes can be used 
to explore different dynamical regimes of the system, such 
as superfiuid, collisional, or collisionless regimes, for both 
Bose and Fermi statistics [1-4]. On the other hand, the 
mode frequencies allow us to probe the equation of state 
(EOS) of the system, including its temperature depen- 
dence. Major benefits result from the high accuracies 
attainable in the measurements of collective frequencies, 
which often enable refined investigations of subtle inter- 
action effects. 

The many-body physics of unitary Fermi gases, i.e. 
two-component Fermi gases with infinite scattering 
length, has attracted tremendous interest over the past 
decade [2-4]. The unitary Fermi gas is characterized by 
strong interaction effects in the EOS [5-8] and reveals a 
unique universal thermodynamic behavior [9]. Further- 
more, at finite temperature, the strong interactions favor 
the collisional hydrodynamic regime, differently from the 
common situation in weakly interacting Bose gases. The 
low-frequency modes of a harmonically trapped Fermi 
gas have been the subject of intensive experimental [10- 
12] and theoretical (see [3] and references therein) ef- 
forts. The temperature dependence has been studied in 
Refs. [13-15]. Remarkably, at unitarity all modes ob- 
served so far turned out to be insensitive to the dif- 
ferent nature of a superfiuid and a classical gas, with 
their frequencies remaining independent of temperature 
throughout the hydrodynamic regime. Previous exper- 



iments have demonstrated the crossover from hydrody- 
namic to collisionless behavior, which typically occurs for 
temperatures approaching the Fermi temperature, with- 
out giving any further information on the regime of lower 
temperatures where the gas is deeply hydrodynamic and 
the superfiuid phase transition occurs. 

In this Letter, we report a joint effort of theory and ex- 
periment on higher-nodal collective modes in the unitary 
Fermi gas. We present a ID hydrodynamic approach to 
describe axial modes in a trapped 'cigar-shaped' cloud. 
Our experimental results confirm the predicted intrin- 
sic sensitivity of higher-nodal modes to the EOS in the 
low-temperature regime, above and even well below the 
superfiuid phase transition. 

The macroscopic dynamic behavior of a superfiuid is 
governed by the Landau two-fluid hydrodynamic equa- 
tions [16] holding in the deep collisional regime lot <C 1 
where r is a typical collisional time and to is the frequency 
of the relevant sound mode in the trap. Landau's equa- 
tions consist of the equation of continuity for the total 
density, the equation for the velocity of the superfiuid 
component, the equation for the entropy density and the 
equation for the current density. The physical ingredients 
entering these equations are the equation of state and the 
superfiuid density. At zero temperature they reduce to 
the irrotational hydrodynamic equations of superfluids, 
above the critical point to the usual hydrodynamic equa- 
tions of normal fluids. Below T c these equations describe 
the propagation of first and second sound, the former 
being basically a density wave, with the normal and su- 
perfiuid components moving in phase, the latter being 
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a temperature or entropy wave. For weakly compress- 
ible fluids the coupling between first and second sound 
is small [16]. This is the case of superfluid helium and 
also of the unitary Fermi gas [17]. Since in the present 
work we investigate the collective oscillations of density 
(first sound), we simplify the search for the solutions of 
Landau's equation by requiring that the velocity fields of 
the normal and superfluid components are equal. Under 
this approximation the equations of motion involve only 
the EOS, the superfluid density playing a role only in the 
propagation of second sound. 

In the presence of a trapping potential, the solution of 
Landau's equations is highly nontrivial, due to the inho- 
mogeneity of the density profile, and thus far has been 
calculated only for the simplest case of isotropic trap- 
ping [17, 18]. Since the experimental excitation and ob- 
servation of these modes are more easily accessible with 
very elongated traps, in the following we focus on such 
configurations. In [19] it was shown that under suitable 
conditions of radial trapping one can derive simplified 
ID hydrodynamic equations starting from the more gen- 
eral 3D Landau's equations. The basic point for such a 
derivation is the assumption that both the velocity field 
v z along the long axis and the temperature fluctuations 
during the propagation of sound do not depend on the 
radial coordinates. This ID-like hydrodynamic formu- 
lation is justified under the condition that the viscosity 
and the thermal conductivity are sufficiently large to en- 
sure, respectively, the absence of radial gradients in the 
velocity v z and the temperature. The condition for the 
viscosity can be recast in the form r\ 3> p n i^ where p n \ is 
the normal ID mass density, obtained by radial integra- 
tion of the 3D normal density. An analogous condition 
holds for the thermal conductivity. These conditions are 
better satisfied in the presence of tight radial confine- 
ment and for the lowest frequency oscillations. One can 
estimate the values of shear viscosity rj using the exper- 
imental data of [20]. For the actual conditions of our 
experiments both sides of the inequality are of the same 
order of magnitude and consequently the full applicabil- 
ity of the ID hydrodynamic formulation can be justified 
only a posteriori. Violation of the ID condition would re- 
sult in a damping of the collective oscillations so that, to 
the extent that the observed damping is small, we expect 
that the assumption of velocity field and temperature be- 
ing independent of the radial coordinate is a reasonable 
ansatz for our variational approach. 

Under the above assumptions and focusing on the uni- 
tary Fermi gas, the equation for the velocity field, char- 
acterizing the density oscillations of the gas in a highly 
elongated harmonic trap, takes the form (see [21] for a 
complete and systematic derivation) 
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tional to e~ ZUJt . Here, uj z represents the trap frequency 
along the axial direction z, and m is the atom's mass. 
Equation (1) explicitly points out the crucial role played 
by the equation of state Pi(ni 7 T), where P\ = j Pdxdy 
is the "ID pressure" (having units of force) and n\ = 
J ndxdy is the atom number per unit length. In order to 
derive Eq.(l) we have explicitly used the adiabatic result 
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7/5Pi holding at unitarity at all tem- 



peratures, where s\ = (1/ni) / sdxdy is the entropy per 
atom with s the entropy density and we have assumed 
the validity of the local density approximation along both 
the axial and radial directions. 

At zero temperature and in the classical limit of high 
temperature the hydrodynamic equation (1) admits ana- 
lytic solutions of polynomial form v z = a\ t z k +a,k-2Z k ~ 2 + 
with integer values of k. At T = 0, where Pi (ni)/m = 
(2/7) [/xo — (^/2) mu} lz 2 ], with /jo the chemical potential 
at the center of the trap, the frequency of the fc-th mode 
is given by 



u 2 = \(k + l)(fc + 5)c 



(2) 



In the classical limit, where P\jn\ = ksT, one finds the 
different k dependence 



uj 2 = i(7fc + 5)cj 2 



(3) 



It is worth noting that Eqs. (2) and (3) coincide for 
k = (center of mass oscillation) and k = 1 (lowest ax- 
ial breathing mode), while they predict different values 
for the higher nodal solutions. One can actually prove 
that not only the frequency of the center of mass but, at 
unitarity, also the frequency of the lowest axial breathing 
mode are independent of temperature, corresponding to 
an exact scaling solution of the two fluid hydrodynamic 
equations [21]. It then follows that only the k > 2 modes 
exhibit a temperature dependence. 

In order to provide a simple quantitative prediction 
for the temperature dependence of the mode frequencies, 
we have developed a variational approach to the solution 
of the hydrodynamic equation (1) with the ansatz v z = 
C12Z 2 + ciq for the k = 2 mode. This ansatz reproduces 
exactly the frequency of the k = 2 mode in the limits of 
T = and high T. For intermediate temperatures we 
obtain [21] 
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where t 2 = MqM^/M 2 . We have introduced the dimen- 
sionless moments 



dx^n -x) {i+1)/2 n(x)X 



(5) 



where we have considered oscillations in time propor- 



of the 3D number density n(x), where x is the local chem- 
ical potential times j3 = 1/ksT and At is the thermal de 



3 



Broglic wavelength. The temperature dependence of the 
moments Me can be evaluated using the experimentally 
determined EOS [8]. Approaching the classical regime, 
for j3(x < —1.5, the virial expansion of the EOS [22] holds 
[8] and allows us to extend the integral to — > — oo. 
At low temperatures, corresponding to /3/x > 4, the EOS 
is governed by phonons and is solely determined by the 
Bertsch paramter £. The error of the quantity above 
resulting from the error in the density EOS is less than 
1%. We have checked that our variational predictions 
are practically indistinguishable from the exact numer- 
ical solution of Eq. (1). Starting from the equation of 
continuity one can also calculate the shape of the density 
oscillations of each mode, proportional to d z {ri\V z ). 

Experimentally, we prepare an ultracold, resonantly in- 
teracting Fermi gas by evaporating a two-component spin 
mixture of fermionic 6 Li in an optical dipole trap [23]. 
The atomic cloud contains typically N/2 = 1.5 x 10 5 
atoms per spin state, and the magnetic field is set to 
834 G, right on top of the well-known broad Fcshbach 
resonance [24]. For the lowest temperatures, the trap- 
ping beam (wavelength 1075 nm) has a waist of 31 urn, 
the trap depth is about 2 fiK, and the axial and ra- 
dial trap frequencies are uj z — 2ir x 22.52(2) Hz and 
u> r = 2ir x 473(2) Hz, respectively. For experiments 
at higher temperatures, the beam waist is increased to 
38 urn, and deeper traps are used (up to 16 /iK depth) 
with trap frequencies of up to uj z = 2tt x 23.31(3) Hz and 
uj r = 2ir x 1226(6) Hz. The corresponding Fermi tem- 
peratures Tp = HfiNuirUJz) 1 ^ 3 /ks vary between about 
0.8 and 1.5 uK. We point out that essentially perfect 
harmonic confinement along the long trap axis (z-axis) 
is ensured by the magnetic trapping that results from 
the curvature of the magnetic field used for Feshbach 
tuning [23]. Also anharmonicities in the radial confine- 
ment remain negligibly small. To probe the ultracold gas 
we record one-dimensional axial density profiles ni(z) by 
near in situ absorption imaging [25]. 

The temperature T of the gas is set by controlled heat- 
ing, always starting from a deeply cooled cloud (T/Tp « 
0.1). In the low-temperature range (T < 0.2 Tp), we 
simply introduce a variable hold time of up to 4 s in 
which residual trap heating slowly increases T. Higher 
temperatures are reached by parametric heating, mod- 
ulating the trap power at about 2uj r , and introducing 
a sufficient hold time to reach thermal equilibrium be- 
tween the different degrees of freedom. We character- 
ize the resulting temperature in a model-independent 
way that docs not require any a priori knowledge of 
the EOS. Based on the virial theorem [26] we intro- 
duce the dimensionless parameter E/E , which repre- 
sents the total energy E = 3mu} z dzz 2 n\(z) normal- 
ized to the energy of a noninteracting, zero-temperature 
Fermi gas, Eq = jNksTp. For a given EOS, the energy 
scale (E/Eq) can be converted to a temperature scale 
(T/Tp). Alternatively, we obtain the cloud's tempera- 



ture by fitting the experimental profiles n\(z) [7, 27] with 
T-dcpendent profiles for a given EOS. For both methods, 
we use the EOS from Ref. [8]. We note that the temper- 
atures obtained by both methods in general show satisfy- 
ing agreement with each other. At very low temperatures 
the latter method shows a trend to give slightly lower val- 
ues of T (up to ~10%), which indicates small systematic 
uncertainties of our measurements. 

We selectively excite axial modes of order k by using 
a resonant excitation scheme. As illustrated in Fig. 1 (a) , 
a repulsive 532-nm laser beam perpendicularly intersects 
the trapping beam, with its position and size chosen in 
a way to provide best mode matching. Typically, the 
excitation pulse contains 8 cycles of sinusoidal modula- 
tion with a half-cycle sine envelope, and the maximum 
potential height of the excitation beam is kept to about 
O.OlfcsTp. The power, length and shape of the exci- 
tation pulse are optimized in order to resonantly drive 
the desired small-amplitude oscillation. The amplitude 
of the corresponding density modulation stays well below 
3% of the central density of the cloud. 

We record axial density profiles n\(z,t) of the excited 
cloud for various time delays t after the excitation pulse. 




-2-1012 

FIG. 1: (Color online) Probing a higher-order first sound 
longitudinal mode at the example of k = 2. In (a), we il- 
lustrate the basic geometry of exciting the optically trapped 
cloud with a weak, power-modulated repulsive laser beam, 
which perpendicularly intersects the trapping beam. In (b), 
the experimental mode profiles (data points) are compared 
to theoretical curves based on the experimental EOS from 
[8] (solid lines) for two different temperatures. The corre- 
sponding cloud profiles in (c) are analyzed to extract the 
temperatures (see text). The solid lines show the density 
profiles obtained from the EOS [8] with T/T F = 0.10 and 
0.45. The parameter Ztf = £ w^ 1 ^/IIibTf /m represents 
the Thomas-Fermi radius of the zero-T interacting gas. 
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We then perform a Fourier transform (FFT). The result- 
ing function h(z,u>) reveals the collective mode spectrum 
with eigenfrequencies and the corresponding spatial 
mode profiles hk(z). We extract the precise frequency of 
a particular mode k by projecting ri\(z, t) onto the spa- 
tial profile hk{z) and analyzing the resulting oscillation 
in the time domain [25]. The high signal-to-noise ratio 
results in very low statistical uncertainties in the permille 
range. 

In this way, we study the longitudinal modes with 
k = 0,1,2. The measured frequency wo of the sloshing 
mode (k = 0) is an accurate measure of the axial trap 
frequency (uJk=o =w z ), and is therefore used for normal- 
ization purposes. The axial compression mode (k = 1) 
has been studied in previous work [11, 28]. Here, in the 
full temperature range explored (0.1 < T/Tp < 0.5) we 
observe its frequency very close to Wfe=i = a/12/5 us z . 
Deviations from this value remain below 0.3% and no sig- 
nificant temperature dependence is observed. This con- 
firms that this mode is insensitive to the temperature as 
long as the gas stays hydrodynamic. For the higher-nodal 
mode with k = 2, we observe the expected T-dependent 
frequency variations. Damping increases as compared to 
the k = 1 mode, but remains sufficiently small to observe 
many oscillations and thus to accurately determine the 
mode frequency. Typical observed damping times [25] 
are 2 s at 0.5 T c , 0.2 s at T c , and 0.12 s at 2T C . We note 
that the k = 3 mode [25] shows very similar behavior, 
with larger frequency variations but faster damping. 

In Figs. 1(b) and (c), we show examples for the spatial 
profiles of the k = 2 mode for two different temperatures 
T/Tp = 0.10 and 0.45 along with the corresponding un- 
perturbed density profiles of the cloud. The comparison 
of the experimental data (data points) with the theoret- 
ical results based on the experimental EOS in [8] (solid 
lines) shows excellent agreement. 

Figure 2 presents the comparison between the exper- 
imental and theoretical frequencies for the k = 2 mode. 
In (a) the normalized mode frequencies ujk=2/w z are plot- 
ted versus the energy parameter E/Eq, while (b) displays 
the same data on a temperature scale T/Tp. The experi- 
mental data confirm the pronounced temperature depen- 
dence of the mode frequencies as predicted by our theory 
based on the EOS of Rcf. [8] (solid line). In comparison, 
the disagreement with the dependence that would result 
from the EOS of the ideal Fermi gas (dashed line) high- 
lights the important role of the EOS. At the lowest tem- 
peratures (T/Tp w 0.1) the frequency lies close to the 
T = superfluid limit, but already shows a significant 
down-shift amounting to almost 1%. At the highest tem- 
peratures (T/Tp w 0.45) our data show a clear trend to 
go below the asymptotic high-temperature value, i.e. the 
classical hydrodynamic limit. The corresponding non- 
monotonic temperature dependence can be understood 
based on the first-order correction to the EOS resulting 
from the virial expansion at high temperatures. 
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FIG. 2: (Color online) Comparison between experimental and 
theoretical first sound frequencies of the k = 2 mode. In (a) 
the experimental data are plotted versus the energy param- 
eter E/Eo, while in (b) we use a temperature scale T/Tp. 
The theoretical curves (solid lines) are based on the EOS of 
Ref. [8]. This EOS is also used to extract T/Tf from the mea- 
sured profile in two different ways: The filled symbols in (b) 
result from a direct conversion of E/Eq to T/Tf, while the 
open symbols result from fitting the cloud profiles (see text). 
For comparison, we also show the mode frequencies (dashed 
curves) that would result from the EOS of the ideal Fermi 
gas. The thin horizontal dashed lines mark the zero-T super- 
fluid limit (lj/lu~ = a/21/5) and the classical hydrodynamic 
limit (uj/uj z = a/19/5) according to Eqs. (2) and (3), respec- 
tively. In (a) the dashed vertical line indicates the T = 
ground state with E/Eo = a/£ = 0.613(3), while the dash- 
dotted vertical lines in (a) and (b) indicate the critical energy 
E c /E = 0.934(39) and temperature T C /T F = 0.223(15). 



In conclusion, our combined theoretical and experi- 
mental work on higher-nodal axial collective modes of a 
unitary Fermi gas reveals a pronounced temperature de- 
pendence below and near the superfluid phase transition. 
The observed temperature dependence is a unique feature 
of higher-nodal modes, not present for any other collec- 
tive mode studied in Fermi gases so far. Our theoretical 
approach is based on a ID two-fluid hydrodynamic model 
describing the frequently used elongated 'cigar-shaped' 
trap geometry. The excellent agreement with the exper- 
imental results provides a stringent test for the validity 
of this ID approach and highlights its potential power 
to accurately describe second sound modes. Moreover, 
our measurements provide an independent confirmation 
of the recently measured EOS of a unitary Fermi gas. 
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